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Abstract 

We present a simple and systematic method to calculate the 
Rindler noise, which is relevant to the analysis of the Unruh effect, 
by using the fluctuation-dissipative theorem. To do this, we calculate 
the dissipative coefficient explicitly from the equations of motion of 
the detector and the field. This method gives not only the correct 
answer but also a hint as to the origin of the apparent statistics in- 
version effect. Moreover, this method is generalized to the Dirac field, 
by using the fermionic fluctuation-dissipation theorem. We can thus 
confirm that the fermionic fluctuation-dissipation theorem is working 
properly. 



1 Introduction 



For a uniformly accelerated observer in the flat spacetime, the ordinary 
Minkowski vacuum, which is defined by using the "frequency" with respect 
to the Minkowski time, can be seen as a thermal bath of particles 0. (The 
temperature of the thermal bath is ha/27r where a is the acceleration of the 
observer.) These particles are considered as the Rindler particles which are 
defined by using the "frequency" with respect to the Rindler time || |], U- 
This effect is called the Unruh effect. 

In order to see this effect, it is convenient to introduce a model of a 
"particle detector" called the DeWitt detector . The DeWitt detector is a 
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point-like object with internal energy levels and moves along the world line 
x(t) where r is the proper time of the detector. The internal energy levels 
of the detector are {E^} and corresponding eigenstates are {\Ei)}. That is, 

H D \E i ) = E i \E i ), (1.1) 

where Hu is the Hamiltonian of the internal structure of the detector. Also, 
the detector has an internal degree of freedom ("monopole moment") Q{r). 
The time evolution of Q{t) in the Heisenberg picture is 

Q(r) = e iHDT/n Q{0) e ~ iHDT/n . (1.2) 

And then, the detector is linearly coupled to a field <f)(x) via this monopole. 
The interaction is described by the interaction Lagrangian 

L^t = g(r)0(r), (1.3) 

where 0(r) = </>(x(r)) is the value of the field along the world line of the 
detector. (To be precise, we need the adiabatic switching.) 

The measurement process by this detector is as follows. Suppose that, 
at r = — oo, the detector is in the ground state \E ) and the field is in the 
Minkowski vacuum |0m)- After the detector- field interaction is switched on, 
the detector would not remain in \Eo), but would make a transition to an 
excited state. It means that the detector "detects" some particles. 

The transition amplitude for the detector-field system to be found in 
| .Ei, ip) at r = oo is given by first order perturbation theory as 

/oo 
drQ{r)(f>{T)\E ,0 M }. (1.4) 

(Here, it has been assumed that the matrix element of Q is sufficiently small 
enough for the perturbation theory to be appropriate.) By Eq. ( |1.2[ ), this can 
be written as 

/oo 
dTe ^-E )r/K m{T)lQM y ^ 
-oo 

Therefore, after summation over all the final states of the field the tran- 
sition rate (the transition probability per unit proper time) from E to E\ 
is 

K^IQWlEo)! 2 ^^^), (1-6) 
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where ^ 

= [ d{r- t') e- iuJ ^ T,) g(r - r') (1.7) 

and 

(7(r-r') = (OM|0(r)0(r')|OM). (1.8) 

Thus, the transition rate of the DeWitt detector is proportional to the 
"response function" J-"(u}) which depends only on the field and the world 
line of the detector but not on the internal structure of the detector. Note 
that we may regard g{r — r') as a kind of noise, "the quantum noise in the 
Minkowski vacuum along the world line x(r)" and the response function as 
the power spectrum of this noise. 

For the Unruh effect, the relevant noise is the Rindler noise, where x(t) is 
a uniformly accelerated world line. In two dimensions, the power spectrum of 
the Rindler noise (the response of the detector) is exactly those of the thermal 
noise in the thermal bath. However, in other dimensions, it was shown by 
Takagi || that there are some differences between the Rindler noise and the 
thermal noise. Specially, the Rindler noise exhibits the phenomenon of the 
apparent inversion of statistics in odd dimensions. 

In this paper, we present a simple and systematic method to reproduce 
these results by using the fluctuation-dissipation theorem |?J §, || which is 
a basis of statistical mechanics for irreversible processes when the systems 
are slightly away from thermal equilibrium. This theorem states the relation 
between the spontaneous fluctuation of fields in thermal equilibrium and the 
irreversible dissipation. Although the fluctuation-dissipation theorem has 
been formulated by various authors, we adopt the formulation by Callen and 
Welton [0], because their formulation is intuitively understandable and ap- 
pealing. They showed that a general form of fluctuation-dissipation theorem 
covers a wide range of phenomena such as the Einstein relation for Brownian 
motion, the Nyquist formula for voltage fluctuation in conductors, and the 
Planck distribution for photons. 

Of course, the "fluctuation-dissipation point of view" was already pointed 
out JTU], |5j . However, the previous discussions concentrated on the vacuum 
expectation value of the commutator and anticommutator of the fields in 
connection with the Huygens' principle. Although these commutators are 
related to the dissipation, the notion of "dissipation" is not quite clear. In 
addition, the calculations involved in obtaining these commutators are quan- 
tum and essentially the same as those to obtain the "fluctuation" directly. 
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In contrast, in our application of the fluctuation-dissipation theorem to 
the Unruh effect, we calculate the dissipative coefficient explicitly from the 
equations of motion of the detector and the field. Then, by virtue of the 
theorem, we can immediately obtain the Rindler noise. All the calculation 
we need is classical and thus our calculation is completely different from the 
previous ones. By this calculation, we can get not only the correct answer 
but also a hint as to the origin of the apparent statistics inversion effect. 

Moreover, our method is generalized to a Dirac field, by using the 
fluctuation-dissipation theorem of a fermionic operator ||11|| . In the context 
of condensed matter physics, the fermionic fluctuation-dissipation theorem 
may not be directly applicable, because we usually measure bosonic quan- 
tities such as voltage or electric current. In this paper, we show that the 
fermionic fluctuation-dissipation theorem is indeed applicable to the Unruh 
effect and works properly. 

To our knowledge, the present approach to the Unruh effect by using the 
bosonic and fermionic fluctuation-dissipation theorem as a corner stone has 
not been discussed before. 

2 Rindler Noise of Real Scalar Field 
2.1 Thermal Noise 

In order to see how the fluctuation-dissipation theorem works not only for 
the thermal noise but also for the Rindler noise to be discussed later, we first 
consider the thermal noise. The system consists of a real scalar field and a 
detector in n dimensional flat spacetime. The detector is at rest and linearly 
coupled with the scalar field through an internal degree of freedom Q(t), 
there. The scalar field is initially in thermal equilibrium at temperature T. 
Thus, the action of the total system is 




(2.1) 



where 
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= J dtQ(t)<f>(x ,t), 



dtdx - 



n-1 



i=l 



(2.3) 
(2.4) 



and Xo is the position of the detector, (x stands for an (n — 1) -dimensional 
vector.) Here, we do not need an explicit form of L(Q,Q) which is the 
Lagrangian for the (unperturbed) detector. 

From this action, we can derive the equations of motion 

+ J dx <j)(x, t) S(x — x ) 



5Q 



9t<P - E + m2( t> - 5 ( f " f o) 



0, 
0, 



where 



(S_So\ 
~ \5Q) 

By Fourier transformations, 



4>(k,u) 

these equations become 

/ SS X 



d fdL\ 
dt [dQj 

dtdx (f)(x, t) e 



'dL_\ 



ikx—iuit 



j dt Q(t) 



-iu>t 



5Q 



+ 



■ikxo 



o, 

0. 



(2.5) 
(2.6) 

(2.7) 

(2.8) 
(2.9) 

(2.10) 
(2.11) 



( — uj 2 + \k\ 2 + m 2 ) (f)(k, u) - Q(uj) e tk *° 

After elimination of <f)(k, u) from these equations, one finds the effective 
equation of motion for the detector 



5Q 



+ K n (u) Q{uj) = 0, 



(2.12) 



where 



K n (u) 



dk 



1 



(2tt 



\n— 1 



2 2 ~ n 7T^" ' x 



— uj 2 + \k\ 2 + m 2 ± it 

K n ~ 2 dK 



r(2fi) Jo -u 2 + k 2 + m 2 ±ie 



(2.13) 
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Note that we have added the ±ze terms to the denominator of the integrand 
to avoid the singularity. The sign is + for uj > and — for u < due to 
causality. (The integrand, as a function of u, must be analytic in the lower- 
half plane of the complex u plane.) Thus, K n (uj) has the imaginary part. 
This gives the friction term. We note that the conventional definition of the 
friction term contains the time derivative of Q(t), which is converted to iu 
by Fourier transformation. So, the conventional definition of the dissipative 
coefficient R n {uS) is given by 

lmK n (u) = -ujR n (u). (2.14) 

The real part of K n (u) would diverge for higher n. But, fortunately, for our 
application of the fluctuation-dissipation theorem, we only need the imag- 
inary part. (The divergence of the real part of K n (oS) would generally be 
renormalized by the potential of the detector.) 
By using 

— !— = P - =F in6(x) (2.15) 
x ± ie x 

and 

S( x 2 - a 2 ) = ^— [6(x + a) + 5{x - a)} , (2.16) 
2\a\ 

one finds 

lmK n (u) = T 2 L^i! (y^-m 2 ) n ~ 3 6(lu 2 - m 2 ). (2.17) 
Thus, the dissipative coefficient R n (ui) is 

Specially, for the massless case, m = 0, this becomes 

finM = ^ W n ~\ (form = 0). (2.19) 

By the fluctuation-dissipation theorem 0, this dissipation suggests the 
fluctuation of the scalar field if in thermal equilibrium. The fluctuation of 
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the scalar field is defined by 

(0(x o ,t)0(xo,*)}/3 = Tr 



-pn 



(f)(x ,t)<f)(xo,t) /Tr 



dt'((l){xQ,t)<f){x ,t'))pS{t-t') 



— / duj F n (u) 

Z7T J-oo 



where 



F n (u) = I ^'e-^ t - t, )(0(f o ,t)0(xo,O)/3 



(2.20) 



(2.21) 



is the power spectrum of the thermal noise. Then, the fluctuation-dissipation 
theorem M says that, 



du F n (u>) 



duo 



dio 



1 hu) 

-hu H — — 

2 - 1 



huRn(\u}\) 



• -x © 
2 2 -7T- 



9 r> 3 ~" 



phui i 

hujluj 



n— 4 



do; 



z,phu) 



(2.22) 

(for m = 0).(2.23) 



This is perfectly consistent with the previous result. (See, for example, 
Ref.i.) 



2.2 Rindler Noise 

Next, we consider the Rindler noise. The system again consists of a real 
scalar field and a detector in n dimensional flat spacetime. However, the 
detector is now uniformly accelerated, i.e., the world line of the detector is 



x°(t) = a^sinhar, (2.24) 

x l (r) = a -1 cosh ar, (2.25) 

x\t) = const., (i = 2,---,n-l) (2.26) 

where r is the proper time of the detector and a is the acceleration. Therefore, 
it is convenient to take the Rindler coordinates |T2J 

x° = £sinhr/, (2.27) 

x 1 = £cosh??, (2.28) 

x i = x\ (i = 2,---,n-l) (2.29) 
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In these coordinates, the world line of the detector is 



ij(t) = ar, (2.30) 
£(r) = a~\ (2.31) 
x\t) = const. (i = 2, ■ -.,n - 1) (2.32) 

The detector is again linearly coupled with the scalar field through an internal 
degree of freedom Q(t), there. However, the scalar field is initially in the 
ground state in this case. The action of the total system is 

S = S (Q) + S int (Q, <P) + S ((/>), (2.33) 

where 

S (Q) = J dr L(Q,Q), (2.34) 
Sint(Q,(p) = J drdrjd^dx Q(r) 4>(£,x,r]) 8(£ - cT 1 ) 5(r) - ar) 5(x - x ) 
di~ Q{j) (j)(a -1 , x , or), (2.35) 



S (4>) = J drjd^dx ~£ 



75(M 2 -(^) a -EW) 2 - m V 

S i=2 



and a^o is the position of the detector, (x stands for an (n — 2)-dimensional 
vector.) 5*o (0) is the same as Eq.( |2.4j ) but written in the Rindler coordinates. 
Again, we do not need an explicit form of L(Q,Q) which is the Lagrangian 
for the (unperturbed) detector. 

The equations of motion derived from this action are 



5_So 
5Q 



+ J drjd^dx 0(£, x, 77) 

x 5(£-a- 1 )5(r ] -aT)5(x-x ) = (2.37) 



and 



- ^Jdr Q(r) 5(i - a' 1 ) 5(t] - ar) 5(x - x ) = 0, (2.38) 
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where 

f6S n \ d (dL\ (dV 



\SQ J dr\dQj \dQ / 
In this case, we consider the transformations, 



(2.39) 



I f 00 dt r 

$(i/,fe,f2) = ~ I -r / d7]dx \r2n7^mhwu K iu (M k i) 

71 Jo £ J 

xe ikx - in n4>(Z,x,r)), (2.40) 
Q(w) = / dr Q(r) e~ iuJT , (2.41) 



where 

M fc = ^m 2 + |fc| 2 (2.42) 
and K iu (z) is a modified Bessel function of imaginary order which satisfies 

\fl& + z Tz~ {z2 ~ u2) ] Kiv{z) = °- (2 ' 43) 

(Note that Q is a frequency with respect to the Rindler time rj and u is 
a frequency with respect to the proper time of the detector r.) By the 
orthogonality relation 

1 f°° dx T ^ . . T ^ . . 5(u — i/) , , . 

- / — K ¥ (x)K iv (x) = ^ J , 2.44 
7r jo x 2z/smh7rz/ 

we can write the inverse transformation of Eq. Q2.40 ) as 

] roc r 

</)(£, x,rj) = / du / dfidfc V 2^ sinh 7rz/ K iv (M k £) 

(27T) n_i 7T JO •/ 

x e- ikx+tQ ^(u,k,n).(2A5) 

Then, the equations of motion become 

fss \ 1 



5Q J (2n) n - 2 na Jo 



r 

dv J dk V2u sinh ttv Ki V {M k /a) 

Xe- ikXo <$>(v,k,uj/a) = 0,(2.46) 



-Q 2 + v 2 ) $(z/, k, Q) V2u sinh nv K iu (M k / 



a) 



71 



xe ikx °Q{Qa) = 0.(2.47) 
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After elimination of $(z/, k, Q) from these equations, one finds the effective 
equation of motion for the detector 



'SSo 
5Q j 



+ )C n (u)Q(u) = 0, 



(2.48) 



where 



2 n ~ 2 ir n a 

2 3-n 



00 r 2z/sinh7rz/ 
av \ dk 



7r ^r(^) a Jo 



dv 



(00 /a) 2 + v 2 ± it 
2z/sinh 7rz/ 



[^(M fc /a)] s 



— (uj/a) 2 + v 2 ±it 
x / dKK n - 3 [KUM K /a)} 2 . (2.49) 



Again, we have added the ±ze term to the denominator of the integrand. 
Specially, for the massless case, m = 0, we can use a closed form of integration 
given by the formula 



in— 5 



dx x n ~ 6 [Kp(x)Y 



Thus, one obtains for Eq. (|2.49[) 



r 



71-2 



i 2 



r(n-2) 

xrQ-i + ^r^-i-A*). (2.50) 



a n - 3 r(^; 



1 



27T2r>-2)A) -(cj/a) 2 + z/ 2 ±ie 
where we have used 



r(§-i + ii/) 



|r(*y)| : 



7T 



y sinh 7T?/ 



(y : real). 



(2.51) 



(2.52) 



(Although k integral is absent for n = 2, this expression is also valid for 
n = 2.) By using, 

2 2z „ / 1' 



r(2z) 



2v^F 



r(*)r * + - 



(2.53) 
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one finds that the dissipative coefficient TZ n (u>) is 
TZJu) 



lmK. n (u) 



—uj 
2 l - n 7i V a r ' 



r( 



n—l • 

2 ■ 



UJ 



r(f-i + ^) 



(form = 0) (2.54) 



and that the ratio of this Rindler case to the thermal case Eq. Q2.19Q is given 
by 



r n (u) = TZ n (uj) I R n (uj) 



UJ 


2-n 


r(f-i + ^) 


a 







By using the elementary formulas 



r{z + l) = zT(z), 



r - + z r 



7T 



COS 71" ,2 



(2.55) 
(2.56) 

(2.57) 
(2.58) 



and Eq. (|2.52|) , one finds, for even 



n. 



. n uj 

r[--i + i- 

2 a 



?? 



-2 + 



n 



-3 + 



a 



= dJuj) 



n-2 



r i 



(2.59) 



and, for odd n, 



, n uj 
T(--l + i- 
2 a 



n 



2 

3) + 



uj 



r 2 

2, 



= d n (uj) 



n-2 /, ,\ -1 
UJ 



. 1 UJ 

r ( - + 

2 a 



(2.60) 



11 



Then, the ratio in Eq.( [2.56 ) can be written as 



r n (u) 



(d n (uj) 

\ tanh(7ru;/a) d n (u) 



(n:even) 
(n:odd) 



Note that 



lim r„ 



that is, 



lim lZ r 



(2.61) 



(2.62) 



(2.63) 



By the fluctuation-dissipation theorem, this dissipation means the fluc- 
tuation of the scalar field if in "thermal" equilibrium. Although the scalar 
field is in the vacuum state in this case, we can view it as in thermal equi- 
librium at temperature T = ha/2ir by the "thermalization theorem". There 
are various versions of the theorem. (For a review, see Ref.||.) 

For example, since the positive Rindler wedge R + (a quarter of the 
Minkowski spacetime, x 1 > \x |) is causally disconnected from the negative 
Rindler wedge R- (another quarter of the Minkowski spacetime, x 1 < —\x°\), 
a uniformly accelerated observer who will be permanently confined within R + 
does not concern with the degrees of the freedom associated with Then, 
by tracing out over these degrees from the Minkowski vacuum, the observer 
gets the thermal density matrix at temperature T — ha/2ir |||. More specif- 

ically mm!, 



(0 M |O«|0 M )=Tr 



e -^ aH R } {+) 



/Tr 



(+) 



(2.64) 



where is an operator only for R + and is the Rindler Hamiltonian 
(the generator of rj translation, i.e., the boost) restricted to there. Note 
that, for the uniformly accelerated observer, the generator of r translation is 
aH^K That is, as far as the uniformly accelerated observer is concerned, the 
expectation value in the Minkowski vacuum of the operator can be seen as an 
ensemble average over the density matrix which is the same form as canonical 
ensemble at temperature T = ha/27r. This effective canonical ensemble is 
true even for interacting fields [|l^, |l7j . 

Equivalently, one can see the thermal character from the periodicity of 
the propagator in the imaginary time [|l^, jl5| . Moreover, it was shown by 



Sewell |19[ in the context of the axiomatic quantum field theory that the 
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Rindler noise satisfies the KMS condition |20| (which is the definition of 
the thermal equilibrium with the systems with infinite numbers of degrees of 
freedom) at temperature T = ha/27i for general interacting field of any spin 
in any dimension. The intuitive explanation of Sewell's theorem is found in 
Ref.0. 

From these various versions of the thermalization theorem, one might 
conclude without calculations that the Rindler noise is the same as the ther- 
mal noise at temperature T = ha/2-K. However, these results only mean that 
the uniformly accelerated observer would see the Minkowski vacuum as the 
thermal equilibrium state but do not mean that the uniformly accelerated 
detector would respond in the same way as it would do at rest in the thermal 
bath, as was emphasized by Unruh-Wald j|] and Takagi ||. The thermaliza- 
tion theorem, itself, says nothing about the noise. In fact, the Rindler noise 
is different from the thermal noise in several points. 

Now, we shall calculate the Rindler noise. By using Eq. fl2.54| ), the ther- 
malization theorem and the fluctuation-dissipation theorem, we can take a 
shortcut to calculate it. The fluctuation of the scalar field is defined by 



(O m |<Mt)0(t)|O m > = — 

Z7T 



du F n {u), 



(2.65) 



where 0(r) = 0(£(r), x(t), T)(t)) is the value of the field along the world line 
of the detector and 



FJu) 



dr' 



j(t-t') 



(O M |0(r)0(r')|O M ) 



(2.66) 



is the power spectrum of the Rindler noise. Then, by using Eq. (|2.64j) , the 
fluctuation-dissipation theorem says that, for the massless case, 



/oo r 
du T n {u) = 4 / 
-oo JO 



oo 

4 / du 





-hu + 

2 

huK n {\u\) 



TZJu) 



2 / dcv 

- 1 

2 / du — — r n (|o;|) 



-oo 
oo 



du 



huR n (\u\) 
ofthu) ^ i Vi 



d n (u), (for m = 0), (2.67) 
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where 



d 2 {u)) 
d 3 (cu) 

d^{ui) 

d 5 (u) 



1 x 



\u\ 



1 + 



l + (- 

00 



a 
a 



x 



10 



\u\ 



(2.68) 
(2.69) 

(2.70) 

(2.71) 

(2.72) 



This result including the statistical inversion in the denominator of Eq. (|2.67|) 
is consistent with that of Takagi ||, except for a minor difference in d n (u) 
for odd n. (We believe that our result is a correct one.) Note that, in Ref.[|], 
the case of the complex scalar field was calculated, but the result is the same 
for the real scalar field as will be explained in the next section. 

We have derived Eq. (|2.67 ) on the basis of Eq.( [2.64j ) and the fluctuation- 
dissipation theorem of a bosonic operator. This clearly shows that the statis- 
tics inversion is an "apparent" one not based on the basic principle of statisti- 
cal mechanics but based on the "temperature" dependence of the dissipative 
coefficient. 



3 Generalization to Complex Scalar Field 

Now, we generalize the analysis to that of a complex scalar field. The sit- 
uation is the same as the real scalar case, except that Q and <fi become 
non-hermitian. 



3.1 Thermal Noise 

For the thermal noise, Eqs. ([2.3|) and 



are replaced by 



Sint{Qy 



S ( 



dtdx 



Q(t)(f>(x,t) + Q\t)(f)\x,t)] S(x-x 



dt 
dtdx 



Q{t)<f>{x ,t) + Q\t) 4t(x ,t) 



n— 1 



|9 t 0| 2 -£|c^| 2 -m 2 |0p 



i=l 



(3.1) 
(3.2) 
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Thus, the effective equation of motion for the detector is 



5Qt 



+ K n (u)Q(Lu) = 0, 



(3.3) 



where K n {u) is the same function as the real scalar case Eq.( 2.13|) . Therefore, 
the dissipative coefficient is the same as the real scalar case Eq. (|2.19|) . 

Then, we can use the fluctuation-dissipation theorem of a non-hermitian 
operator [ITT]. The fluctuation of the field is now defined by 



(<j>(£o,t)$(x Q ,t))p 



1 r°° 



2tt 



du F n (u), 



where 

F n (u) = J dt' e-^-^^x^^t'))^ 
and the theorem says that 



(3.4) 
(3.5) 



/oo poo fujj , 

du F n (uo) = 4 du — 
-oo JO Z ^ 



1 + 



1 



1 



+ 



00 , TiujRJ-u) 
du 



2 2 - n 7i — 



Tiu)\ui 



n— 4 



oo eP hu} - 1 



-Rn(^) 

(3.6) 

(form = 0) (3.7) 



3.2 Rindler Noise 

Similarly, for the Rindler noise, Eqs. (|2.35|) and ( |2.36|) are replaced by 



S int (Q,<j)) 



drdr]d£dx Q{r) <fi(£, x, rj) + Q J (r) T (£, x, rj) 

x 8(£ — cT 1 ) 5(r] — or) 5(x — x ) 



dr 



Q(r)(f)(a 1 , x , or) + Q\t) $ (a 1 ,x ,ar) , (3.8) 



drjd^ dx £ 



n-l 



2| / (2 



i=2 



(3.9) 
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and the effective equation of motion becomes 



5Qt 



K. n {u) Q(u) = 0, 



(3.10) 



where )C n (u) is the same function as the real scalar case Eq. (p.49D . So, the 
dissipative coefficient is the same as the real scalar case Eq. (|2.54 ), again. 
The fluctuation of the field is now defined by 



(O M |0(r)0 t (r)|O M ) = — / du F n (u), 

Z7T J-oo 



(3.11) 



where 



F n {u) = J dr' e-- ^')(O M |0(r)0t( r ')|o M ), (3.12) 
and, by using the fluctuation-dissipation theorem of a non-hermitian opera- 



tor |l Ij and Eq.( |2.64|) , one finds that 

/oo rco YlijJ r 

du T n {u) = 4 / du —\ 
-oo JO Z ^ 



1 + 



1 



du 



-oo 
oo 



du 



+ 

huTZ n (—u) 
huR„(—u) 



1 



-oo - (-l) n 

in agreement with the previous analysis f|. 



K n (u) 

n n {-u)} 

(3.13) 

d n (u). (for to = 0) (3.14) 



4 Generalization to Dirac Field 

Next, we generalize to a Dirac field, using the fluctuation-dissipation theorem 



of a fermionic operator [|TTJ. The situation is the same as the scalar case. 
However, the internal degree of freedom of the detector is now a spinor G 
and the field is a Dirac field ip. 
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4.1 Thermal Noise 

For the thermal noise, the action of the total system is 

s = s (s) + s mt (e,ip) + s (^), 



J dtL(e,e), 

I 



dtdx 



9(t)ij)(x,t)+il)(x,t)Q(t) S(x-x ) 
J dt [e(t) ?p(x , t) + $(x 0: t) G(t)] , 
dtdx %[) [i7 M <9 M — m] ip, 



where 

S (Q) 

s int (e^) 

and 

^ = ^ f 7 , e = ey. 

Here, 7^ are the 7 matrices which satisfy 

The equations of motion become 

/ SlSq 



5Q 



J dx i>(x,t) 5{x — x ) = 0, 



dty -^2^ + m V - l^ 9 ^ + m \ e (*) - x ) = 0, 



(4.1) 

(4.2) 

(4.3) 
(4.4) 

(4.5) 
(4.6) 

(4.7) 
(4.8) 



where "L" denotes the left-derivative. After Fourier transformation and elim- 
ination of the field, one finds that the effective equation of motion for the 
detector is 



5Q 



+ % n He( w ) = o, 



where 



r dk 
J (27T)™- 1 



— CJ7 + k ■ 7 + m 
-uj 2 + \k\ 2 + m 2 ± it 



(4.9) 



(4.10) 
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(The subscript "1/2" denotes the spin of field.) Since the k ■ 7 term is odd 
in k, we can relate Ky 2>n (u) to that of the scalar case Eg. (|2.13|) , 



— uwy + m 



KJlo) 



Ky 2 , n {io) 

Therefore, we define the dissipative coefficient as 



(4.11) 



R 



l/2,n 



-UJ 



-0*7 + m 



Rnipj) 



(4.12) 



where R n {^) is that of the scalar case Eq. (|2.18| ). This extra — sign is a 
convention; it was chosen so that the dissipative coefficient is "positive" in 
the sense that 



Tr 7°-Ri/2,n(^) oc ojR„(oj) > for 10 > 0, 
and this convention is the same as in Ref . [|i~T|j . We also define 



R 



l/2,n\ 



-id 



-R\/2,r 



-0J 



uj'y + m 



Rn(-Uj) 



and 



Tr 7°-Ri/2,n(~ ^) oc toRn{—oj) > for uj > 0. 



(4.13) 



(4.14) 



(4.15) 



Then, we can use the fluctuation-dissipation theorem of a fermionic op- 
erator |TTJ. The thermal fluctuation of the Dirac field is defined by 



{il>(x ,t)il>(xo,t))i3 



Tr 



ip(x ,t)ip(x ,t) /Tr 



-/3H 



dt'(7p(x ,t)t/;(x ,t'))p5(t-t') 
2tt J 



X (lp(x ,t)lp(xo,t'))/3. 

The power spectrum of the Dirac field is defined by 



F 1/2in (u) = A-'Tr 



7 / dt' e-^'^fo,^^,^ 



(4.16) 



(4.17) 
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where A n is the dimension of the 7 matrices. Therefore, 



Tr 



7 O (-0(x o , t)i(}(x , t))p 



2tt 



du Ft 



Then, the theorem [] 1] says that 



/oo 
duj Fi/ 2tn (u) 
-00 



l/2,n 



1 



(4.18) 



e /3hu> _|_ 1 



l/2,n 



+ 



££ <fc '^T ,ft ^(-)] 







-Rl/2,n(— <^)} 



l n J -00 

J -00 + 1 



(4.19) 



This result is consistent with the previous result. (See, for example, Ref. ||.) 
The differences between this case and the scalar case are the power of to 
in the numerator, which comes from the fact that the Dirac equation is 
the first order with respect to the time derivative, and the Fermi distribution 
factor, which is provided by the fluctuation-dissipation theorem of a fermionic 
operator. 



4.2 Rindler Noise 

For the Rindler noise, one might define the Rindler noise of Dirac field as 

(0m|^W(t')|0m). (4.20) 

However, this noise is not stationary, i.e., is not a function of r — t' . Instead, 
we define the Rindler noise of Dirac field as [TJ| 



<0 m |V>(tMt')|0m>, 



(4.21) 



where 



V>(r) 



exp 



1 1 
-aT7 7 



= cosh(ar/2) — 7°7 1 sinh(ar/2) ip{r) 

EE SW(r). 



(4.22) 
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This transformation is the Lorentz transformation, the boost, from the lab- 
oratory frame to the instantaneously comoving frame of the uniformly accel- 
erated observer at r. 

Thus, we write the action in terms of 



0(r) = 5 T e(r). 
Then, the action of the total system is 

s = s (e) + s int (ej) + s {$), 

where 



S (Q) 



drL{e,Q), 
drdrjd^ dx 



(4.23) 
(4.24) 

(4.25) 

(4.26) 



x 6(£ - a~ l ) 5(r] - ar) 5(x - x ), (4.27) 



dr)d£dx ^[i7°7^ + f ^7 + 9 t 



n-l 



(4.28) 



i=2 



This So(ip) is the same as Eq.( |4.4| ) but written in the Rindler coordinates and 
ip. (If we introduce the vielbein, calculate the spin connection in the Rindler 
space and define the spinor with respect to the local Lorentz, we would get 
this action.) 

The equations of motion become 



, 5Q , 



+ / dr/d^dx ip(£, x,r}) 

x 5(£ — a -1 ) 5(rj — ar) 5(x — x ) 



(4.29) 



and 
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2f 



771 



x - fdr Q(r)5^-a- 1 )5(r]~aT)5(x~x ) = 0(4.30) 

In order to deal with the 7°7 1 term of the second equation, we define the 
projection operators 

/ 1±7°7 1 \ 



7± = 



(4.31) 



which satisfy 

7± = 7±> 7+7- = 0, 7+ + 7- 



7dz(7°7 lx 



(4.32) 



By using these projection operators, we can divide tp into two parts 
where ip± = 7±^>. Then, the second equation becomes 

1 / ^±i) 2 ^ ± -^(^±)-E^± + ^± 



e v"~2 

-7± 



X - dr 6(r) <5(£ — a 1 ) 5(i] — ar) <5(a; — x ) 
In this case, instead of K iv {M^), we consider K iv± i (M^) as 



(4.33) 



0. (4.34) 



71 JO £ 



drjdx K iv± i(M k £) 



xe ,b - ffl, l±fc^). (4.35) 



(See, for example, Ref.^TJ.) By the orthogonality relation 



8{n-v) 



JO X 



(=f2iv — 1) cosh txv ' 



(4.36) 
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which can be obtained from Eg. (|2.44f) by the shift iv — > iv± ~, we can invert 
Eg. flOp as 



(2n) n - 1 n 



x 



dvdVLdk (=p2iz/ — 1) cosh7rz/ 
K lu± ,(M k O e- tkx+m ^ ± (u,k,Q). (4.37) 



If we specialize to the massless case, then after a straightforward calcu- 
lation, one finds that the effective equation of motion for the detector is 



, 5Q , 



+ £i/ 2 ,n(w) 9(w) = 0, 



(4.38) 



where 

£l/2,n(^) 



^vrf +1 r(n - 21 



— cj/a + v + ze u/a + v — i 



x 



(2i> — i) cosli7rz/ / n — 3 



x r 



2uj/a — i 
n + 1 



-i-K-y+T 1 ) 

a/ 



— cu/a + i/ + «e tu/a + z/ + z 
(2z/ + z) cosh7rz/ /n — 3 .aA 



x r 



2uo/a + i 
n — 1 



i + iz/ r 



2 a/ 
n + 1 



(7° + 7 1 ) 



1 — iv 



(4.39) 



Note that the singularity on the real axis of v is only at v = u/a. To 
avoid this singularity, we add the +ie term. (The sign of the ie term can be 
determined by the causality, as above. Then, the sign is + for both uo > 
and io < 0.) From this +ie term, we can obtain the dissipative coefficient 



Im/Ci/ 2 ,nM 



— C<J 

2 1 - n vr^ a n ~ l 



r(2f i ; 



coth(7ro;/a) 



r(2±i 



7 



a n 7?.„ + i(a;) coth(7rco>/a) 7 , 



(4.40) 
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where lZ n+ i(u) is that of the scalar case Eq.( |2.54D in n + 1 dimensions, and 



a, 



r(*=i; 



(4.41) 



Here, we have again included the extra — sign as a convention in the definition 
of Eq.flOg). We also define 

^i/2,n(-^) = -Tii/2,n{-u) = a n TZ n+1 (-u) coth(7ru;/a) 7°. (4.42) 

Then, we can use the fluctuation-dissipation theorem of a fermionic op- 



erator [11]. The power spectrum of the Dirac field is defined by 



F 1/2 , n (u) = A^Tr 



7° / dt' e*> ^ (0 M |^(r)^(r')|0 M ) 



and then, 



Tr 



7 (0m|^(t)^(t)|0m) 



A f°° 

ZTT J-oo 



(4.43) 



(4.44) 



Thus, by using Eq. fl2.64 ), the theorem [11] says that 



/°° 4 n f 

du Ty^niu) = ^-Tr 7 / 
-00 l\ n Jo 



2 I 



1 - 



e f3hu> _|_ 1 



1 



e /3hui + 1 



7^1/2,n(-^)} 



2 f°° ftu 

du — -Tr[7°7li/2,n(-^)] 



A, 
2a 



-00 
00 



du 



gfihio _|_ 
huTl n+ i(-u) 



(4.45) 



That is, the power spectrum of the Rindler noise of massless Dirac field in n 
dimensions is proportional to that of massless scalar field Eq. (|2.67 ) in n + 1 
dimensions. From the analysis of the scalar case, we can thus show the 
statistical inversion. This result nicely agrees with that of Takagi ||. 
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5 Conclusion 



We have presented a simple and systematic method to evaluate the Rindler 
noise. We have first calculated the dissipative coefficient explicitly from 
the equations of motion of the detector and the field. Then, by using the 
fluctuation-dissipative theorem, we have obtained the Rindler noise which is 
relevant to the analysis of the Unruh effect. This method is generalized to 



the Dirac field, by using the fermionic fluctuation-dissipation theorem [\ I 



These results are perfectly consistent with previously known results including 
apparent statistical inversion in odd dimensions. 

Although the Rindler noise can be calculated directly from the action of 
the field as Takagi did, we emphasize that there are several advantages in 
calculating the noise indirectly. That is, to introduce the detector, calculate 
the dissipative coefficient and then use the fluctuation-dissipation theorem, 
as we have done in this paper. 

First of all, it contains only the classical calculations expect for the 
fluctuation-dissipation theorem. To calculate the dissipative coefficient, we 
only have to eliminate the degree of freedom of the field from classical equa- 
tions of motion of the detector and field. Thus, the calculations become much 
simpler. 

Next, this method gives a hint as to the origin of the apparent statistics in- 
version effect. While the dissipative coefficient for the thermal noise does not 
depend on the temperature in an ideal setting, those for the Rindler noise in- 
evitably depends on the "temperature" (the acceleration) . This temperature- 
dependence of the dissipative coefficient destroys the simple Bose (or Fermi) 
distribution provided by the fluctuation-dissipation theorem and causes the 
phenomenon of the "apparent" inversion of statistics in odd dimensions [0]. 
Why does this difference happen? From our calculation, the answer is ob- 
vious. Note that all the effect of the field on the detector is represented by 
the influence functional || which is constructed from the action of the field, 
the detector-field interaction and the initial condition of the field. (The in- 
fluence functional for a scalar field in the Minkowski vacuum, coupled to a 
uniformly accelerating DeWitt detector was derived by Anglin f22|.) In the 
thermal case, the temperature is the initial condition. On the other hand, 
in the Rindler case, the temperature is in the detector-field interaction as 
the acceleration. Because it is only the action of the field and the detector- 
field interaction that are needed to obtain the effective equation of motion 
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and thus the dissipative coefficient, the dissipative coefficient for the ther- 
mal noise does not depend on the temperature in our idealized treatment 
but those for the Rindler noise inevitably depends on the temperature. (Of 
course, this point of view would not explain everything about the apparent 
statistics inversion effect.) 

Finally, we can see that the fermionic fluctuation-dissipation theorem [jTT] 



works properly. In the context of condensed matter physics, it is generally 
difficult to use the fermionic fluctuation-dissipation theorem because we usu- 
ally measure bosonic quantities such as voltage or electric current. However, 
we have applied this fermionic version of the theorem to the Unruh effect 
and have obtained the right result including apparent statistical inversion 
in odd dimensions. Thus, we have confirmed that the fermionic fluctuation- 
dissipation theorem is indeed working properly. 

Since the relation between Rindler and Minkowski coordinates is very 
similar to the relation between Schwarzschild and Kruskal coordinates, there 
is a hope that the above method could be applicable to the Hawking radia- 
tion [^, Q also. This would be a future work. 
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